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ABSTRACT 

We investigate where in protoplanetary disks magnetorotational instability operates, 
which can cause angular momentum transport in the disks. We investigate the spatial 
distribution of various charged particles and the unstable regions for a variety of 
models for protoplanetary disks taking into account the recombination of ions and 
electrons at grain surfaces, which is an important process in most parts of the disks. 
We find that for all the models there is an inner region which is magnetorotationally 
stable due to ohmic dissipation. This must make the accretion onto the central star 
non-steady. For the model of the minimum-mass solar nebula, the critical radius, 
inside of which the disk is stable, is about 20 AU, and the mass accretion rate just 
outside the critical radius is 10 -7 - 10 -6 M@ yr" 1 . The stable region is smaller in a 
disk of lower column density. Dust grains in protoplanetary disks may grow by mutual 
sticking and may sediment toward the midplane of the disks. We find that the stable 
region shrinks as the grain size increases or the sedimentation proceeds. Therefore in 
the late evolutionary stages, protoplanetary disks can be magnetorotationally unstable 
even in the inner regions. 
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1. 



INTRODUCTION 



Various infrared and radio observations have revealed the existence of accretion disks around 
T Tauri stars, so-called protoplanetary disks, where planet formation may be going on (e.g., 
Rucinski 1985; Adams, Lada, & Shu 1987). Photometric observations at optically thin millimeter 
wavelengths give disk masses of ~ 0.1 - 0.01 M & (Beckwith et al. 1990), and fitting of spectral 
energy distributions (SEDs) at infrared wavelengths gives disk sizes of tens of AU or 100 AU 
(Adams, Lada, & Shu 1988). Some authors tried to reproduce the SEDs of T Tauri stars 
with models of viscous accretion disks (e.g., Bertout, Basri, & Bouvier 1988), and succeeded 
in deriving the slopes of SEDs at infrared wavelengths for some T Tauri stars. Optical and 
ultraviolet veiling of absorption lines has been discovered for some T Tauri stars, which indicates 
the existence of boundary layers at the interface between a Keplerian accretion disk and a slowly 
spinning star (Lynden-Bell & Pringle 1974). All these suggest the existence of accreting matter 
in protoplanetary disks at least in their early active stages. However, the mechanisms of angular 
momentum transport, necessary for accretion, have not yet been well understood, any more than 
in other astrophysical accretion disks. 

The magnetorotational instability must be the most promising source of anomalous viscosity 
in accretion disks (Balbus & Hawley 1998, and references therein). Balbus & Hawley (1991) 
showed that astrophysical accretion disks are unstable to axisymmetric disturbances in the 
presence of a weak magnetic field. The instability sets in under a broad range of conditions 
and is capable of initiating and sustaining MHD turbulence in accretion disks, as long as the 
approximation of the ideal MHD holds. However, protoplanetary disks are so cold and dense 
that the ionization fraction is very low (Umebayashi & Nakano 1988). Because the coupling of 
some charged particles with magnetic fields is not strong enough, both the processes of ambipolar 
diffusion and ohmic dissipation have to be considered (Nakano 1984). These processes have an 
effect of suppressing the magnetorotational instability. It is necessary to clarify which parts of the 
disks are magnetorotationally unstable in relation to the evolution of protoplanetary disks. 

Blaes &: Balbus (1994) examined the effect of the ambipolar diffusion on the instability, and 
showed that differentially rotating disks are unstable when the collision frequency of neutral 
particles with ions exceeds the local epicyclic frequency. Sano & Miyama (1999, hereafter Paper 
I) investigated the stability of weakly ionized disks including the effect of ohmic dissipation. This 
process dominates the ambipolar diffusion in the regions of higher density, or in the inner parts 
(within about 100 AU) of the disks. 

In paper I, we derived the conditions for the existence of unstable modes by the global linear 
analysis taking into account the effect of the vertical structure of the disks. The wavelength of the 
most unstable mode expected from the local linear analysis is approximately given by 




(1) 



where va = B/(4irp) 1 / 2 is the Alfven speed, f2 is the angular rotational velocity of the disk, and r/ 
is the magnetic diffusivity (Jin 1996; Paper I). If the magnetic Reynolds number R m = v A /riQ is 
larger than unity, the ohmic dissipation is not efficient and the wavelength of the most unstable 
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mode is approximately given by 2-kva/Q, which is the result of the ideal MHD. When R m ^ 1, 
on the other hand, perturbations with shorter wavelengths are stabilized due to the dissipation. 
Thus, as the magnetic diffusivity increases, 2ttt]/va increases and finally becomes the wavelength 
of the most unstable mode. 

The local analysis gives the wavelength Ai oca i of the most unstable mode as a function of the 
local values of va and r/. Thus Ai oca i takes different values at different points in the disk. Our 
results of the global analysis show that the vertical distribution of the wavelength Ai oca i(z) is quite 
important for the understanding of the global stability. In paper I, we have found that the layer 
at a height z is globally unstable when Xi oca i(z) obtained by the local analysis is shorter than the 
scale height of the disk H, or 

Aiocai {*) ~ max 

The condition given by equation (|2|) agrees with the idea that the stability criterion is 
approximately given by the requirement that the shortest wavelength for local instability is smaller 
than the scale height of the disk. 

Thus the ratio \\ oca \/H is a good indicator for the existence of global unstable modes. This 
ratio is sensitive to the magnetic diffusivity, which is determined by the physical quantities in 
protoplanetary disks such as the density, the temperature, and the ionization fraction. Umebayashi 
&; Nakano (1988) examined the radial distribution of magnetic diffusivity on the midplane for 
the so-called minimum-mass solar nebula proposed by Hayashi, Nakazawa, & Nakagawa (1985). 
However, to make use of equation @) as the stability criterion of protoplanetary disks, we 
have to know the vertical distribution of magnetic diffusivity. In this paper we investigate the 
two-dimensional distribution of magnetic diffusivity for various disk models, and reveal where in 
protoplanetary disks the magnetorotational instability is operative. 

Because the recombination processes are more effective in denser regions, the ionization 
fraction is lower at higher density. Thus the inner dense region of the disks would be decoupled 
from magnetic fields except for the innermost region where the thermal ionization is efficient. 
However, even in the inner dense region, the gas near the surface layer of the disks may be well 
ionized because the density is low. Gammie (1996) proposed a layered accretion model which 
has the decoupled region (dead zone) sandwiched by the ionized active layers. The ionization by 
X-rays from the central object would also make a similar structure (Glassgold, Najita, & Igea 
1997). In such situations, the column densities of the active layers are one of the most important 
quantities for the dynamical evolution of protoplanetary disks. We shall calculate the column 
densities of these layers and compare them with the previous works. 

The recombination of ions and electrons on grain surface is the dominant process in the dense 
region. Therefore the characters of dust grains affect significantly the ionization fraction and the 
column density of the active layers. In the previous works grain surface reactions were not included 
because they were concerned mainly in the late phases of the disks when dust grains of the disks 
have already settled toward the midplane. In this paper we investigate the effect of dust grains 
on the abundances of charged particles and the contribution of charged grains to the interaction 
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of the gas with magnetic fields in protoplanetary disks, and show where the magnetorotational 
instability can be operative. 

The plan of this paper is as follows. In §2 we describe the disk models adopted in this 
paper and the magnetic diffusivity of the gas composed of various kinds of charged particles. 
We also describe the reaction scheme which determines the abundance of charged particles in 
ionization-recombination equilibrium. In §3 we show numerical results on the spatial distribution 
of charged particles, the time scale of magnetic field dissipation, and the unstable regions in the 
disks. Because charged grains can be the dominant charged particles at least in some parts of the 
protoplanetary disks, their size and abundance are highly important in determining the magnetic 
diffusivity. Dust grains may grow by mutual sticking. They will sediment onto the midplane if the 
disk is quiescent. The coupling between the gas and magnetic field varies as the disk evolves in 
such ways. In §3 we also examine the unstable regions in some model disks which must correspond 
to different evolutionary stages. In §4 we investigate the column densities of active layers and 
discuss the evolutionary scenario of protoplanetary disks taking account of the magnetorotational 
instability. We also discuss the validity of our assumptions and some related problems in §4. 
Finally we summarize our results in §5. 



2. MODEL AND FORMULATION 

2.1. Structure of Protoplanetary Disks 

Although various models have been proposed for protoplanetary disks, there seem to be no 
really acceptable models yet. In view of such a situation we will investigate on a variety of models 
and find out some firm results. The surface density S(r) of the protoplanetary disks is one of the 
most important quantities in constructing the disk models. We assume a power-law distribution 

where r is the distance from the central star. Hayashi (1981) found that the surface density of the 
primitive solar nebula, which is restored from the mass distribution in the present solar system, 
is well represented by equation (||) with q = 3/2 and So ~ 1.7 x 10 3 g cm' 2 . This model is thus 
called the minimum- mass solar nebula. We examine mainly for the models with q = 3/2 and 
So = 1.7 x 10 3 /s g cm~ 2 in this paper, where /s is an arbitrary parameter. Equation (||) in this 
case must break down at large r because the total mass of the disk should not diverge. In this 
work we examine the regions of r 100 AU. Because the power index q could take different values 
(e.g., Cassen 1994), we will also investigate for the other values of q. 

For the temperature of the disks we also assume a power-law distribution, 

Except for a very early accretion phase, the temperature of the disk is determined by the energy 
balance between absorption of stellar radiation and emission of thermal radiation. After dust 
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grains in the disk have sunk toward the midplane to some extent, most part of the disk is nearly 
transparent to the stellar light and the thermal radiation. The equilibrium temperature in such 
a situation is given by p = 1/2 and Tq = 280 K for the luminosity of the central star L* = 1 Lq. 
We adopt p = 1/2 and To = 280 K even for the opaque regions near the midplane because it may 
be regarded as being in a thermal bath of this temperature. For simplicity we also adopt this 
temperature distribution even for the early opaque stage before significant sedimentation of dust 
grains occurs because the results are not very sensitive to the temperature as long as it is below 



several hundred Kelvin (see §4.5). With p = 1/2 and Tq = 280 K, the sound velocity is given by 



cJr) = = 9.9 x 10 4 — — -f— cm sec -1 , (5 

where k is the Boltzmann constant, p is the mean molecular weight of the gas, and toh is the mass 
of a hydrogen atom. 

In a gravitationally stable disk, the gravity perpendicular to the disk is mainly contributed 
by the central star. The balance of the stellar gravity and the pressure force in the z-direction 
determines the gas density at height z from the midplane as 

p(r, z) = p c {r) exp 



H 2 (r) 



(6) 



where p c {r) is the density at the midplane and 



H(r) = ^ (7) 

is the scale height of the disk. For the Keplerian angular velocity £1 = (GM i /r 3 ) 1 '' 2 , where G is 
the gravitational constant and M* is the mass of the central star, the density at the midplane is 
given by 

/ r \ -11/4 / M n 1/2 / x 1/2 
O'M^xlO-V^) (-1) (£) (8) 

for E given by equation (||) with q = 3/2 and for c s given by equation <^j. 

The local stability analysis shows that a disk rotating around a star with the Keplerian 
velocity is gravitationally stable only when S < Qc s /itG (Goldreich & Lynden-Bell 1965). The 
disk given by equations (||) and (|5|) with q = 3/2 is gravitationally stable when the inequality, 

/ r \ - 1 / 4 /MX 1 / 2 / u \ - 1 / 2 
fc<17.5(— -T71 ) (T7 1 ) Uh) , 0) 



,100 au; \M Q J \2.34 

is fulfilled. Since we are concerned with gravitationally stable disks of radius r <J 100 AU, it is 
sufficient to consider the cases of /s ^ 10. 



2.2. Magnetic Diffusivity 

The magnetic diffusivity 77 is given in terms of the electrical conductivity a c and the light 
speed c by 

"=■£■■ (10> 
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The electrical conductivity is given by 

lv) 2 T v n v 



(7 C 



\^Hv) iy< b v {\X) 



where r v is the viscous damping time of motion relative to neutral particles of a charged particle 
v with mass m„ and electric charge eq u , and n u is the number density of particle v. 

Because the mass of an ion, mi, is generally much greater than the mass of a neutral molecule, 
m n , we have for an ion 

rr 1 = n n (av)i— , (12) 
rri[ 

where n n is the number density of neutrals and {av)\ is the rate coefficient for the collision between 
ions and neutrals averaged over the distribution of the relative velocity. The momentum-transfer 
cross section of a singly charged ion due to electrical polarization force is given by (Osterbrock 
1961) 

2 \ 1/2 
2 



a = 2.41vr ^ I , (13) 



mr 

where a is the polar izability of a neutral molecule, e is the elementary electric charge, m is the 
reduced mass, and v is the relative velocity. Almost all ions at such low temperatures are singly 
charged. We have a = 7.9 x 10 -25 cm 3 for an H2 molecule and a = 2.1 x 10 -25 cm 3 for a He atom. 
By taking the mass of an ion m\ ~ 34 mn, for example, we have (av)i^n 2 ~ 1.8 x 10~ 9 cm 3 s _1 , 
which is independent of v as long as a is larger than the geometrical cross section. Such a rate 
coefficient due to polarization force is called the Langevin's rate coefficient. 

For an electron we have 

Tc 1 = n n (av) e . (14) 

In the previous works, the Langevin's rate was also used for the rate coefficient for the collision 
between electrons and neutrals. However, the collision cross sections obtained by experiments at 
low electron energies are much smaller than a given by equation (13) and is nearly equal to a 
geometrical cross section. In this paper, we use the cross sections recommended by M. Hayashi 
(1981) who compiled the experimental results (see Appendix). 

We neglect the contribution of electron-ion collisions to the electrical conductivity in 
equation (|ll|). The contribution of electron-ion collision is r /r c _i times the contribution of 
electron- neutral collision, where r e _i is the collision time of an electron with ions, and is inefficient 
unless n e /n Q > 1CP 8 and 1CP 4 at T = 10 and 10 3 K, respectively. 

For a grain of any charge, we have 

TpT 1 = n n (av) g — , (15) 

where m g is the mass of a grain. Because the mean free path of neutrals is much greater than the 
size of a grain, we have for a spherical grain of radius a, 

where Cth = (S^T/V^mii) 1 / 2 is the thermal velocity of a neutral molecule. 



-7- 



2.3. Reaction Scheme 

We investigate numerically the densities of various kinds of charged particles in order to 
calculate the magnetic diffusivity 77 given by equation ([ji]). The rate equation for a constituent Xj 
is given by 

= J2 nA^i) + E /%fc"(X>(X fe ) , (17) 

where n(Xj) is the number density of Xj, 7^ is the rate coefficient for forming Xj by ionization 
of Xj, and /3jjfc is the rate coefficient of a two-body reaction including reactions on dust grains. 
Because the gas in protoplanetary disks can be regarded as being in ionization-recombination 
equilibrium (see §4.6 below), equation (|i~7|) is reduced to 

V ^x(Xj) + PijkxQtjWXk) = , (18) 

where x(Xj) = n(Xj)/nn- We solve equation (|T^ ) at each position of the disk to obtain the spatial 
distribution of magnetic diffusivity. 

We consider the regions of the disk where the temperature is less than several hundred 
Kelvin, and thus the thermal ionization is inefficient compared with the ionization by non-thermal 
particles (Umebayashi 1983; Nakano &: Umebayashi 1988). In such regions ions and electrons are 
produced mainly through ionization by cosmic rays and radioactivity, and vanish through various 
recombination processes. Since these situations are quite similar to those in very dense interstellar 
clouds, we adopt a simplified reaction scheme for dense clouds investigated by Umebayashi & 
Nakano (1980, 1990) and Nishi, Nakano, & Umebayashi (1991) with some modifications. Among 
the modifications are updating the reaction rates and adding some species of charged particles 
such as H^" and grains of charge ±3e. We also take into account the electrical polarization of 
grains for collisions with ions and electrons. In the following we describe the reaction scheme 
adopted in this paper. 



2.3.1. Reactions in Gas Phase 

The elements to be considered in our reaction scheme are H, He, C, O, and refractory heavy 
elements. We adopt the solar system abundances for these elements (Anders & Grevesse 1989) 
which are shown in Table [|. Most of the heavy elements have condensed in dust grains and only 
a small fraction of them remain in the gas phase. We denote this fraction as 5\ for the volatile 
elements C and O, and 82 for the refractory heavy elements. We assume for simplicity that all 
grains have the same size unless otherwise stated. Then the grain abundance n g /nH by number 
relative to hydrogen is determined from the total amounts of those elements in grains. 

We can regard that each element in gas phase has been transformed into its molecular form 
because we consider low temperature, high density regions. We assume that all C in the gas phase 
is in CO molecule. Oxygen in the gas phase other than the constituent of CO is in the form of O2, 
H2O, OH, or O. Since the reaction rate coefficients of molecules O2, H2O, and OH with ions are of 
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the same order and the recombination rate coefficients of the resultant ions with electrons are also 
of the same order, we regard the reactions of the most abundant species O2 as the representatives 
of their reactions. We introduce a parameter /q 2 which represents the fraction of oxygen in the 
gas phase in the molecular form O2, H2O, and OH. Because the metal atoms such as Na, Mg, and 
Al have nearly the same rate coefficients for charge-transfer reaction with molecular ions and the 
resultant metal ions also have similar rate coefficients for radiative recombination with electrons, 
we neglect the differences among them and denote them as M collectively. 

We consider the following species of charged particles in our reaction scheme: electron e, 
H + , He + , C + , , Hg~, molecular ion m + (except and H^~), metal ion M + , and charged 
grains. Dust grains in the disks can be charged through the collision with electrons and ions in 
the gas phase. We consider the grains of charge (neutral), ±le, ±2e, and ±3e. Table || shows 
the reactions in gas phase and their rate coefficients (Millar, Farquhar, & Willacy 1997). The 
important reactions of molecular ions m + are dissociative recombination with electrons and charge 
transfer with metal atoms. Because the rate coefficients of these reactions are similar among 
various molecular ions, we regard all m + as a single species. For the rate coefficients of m + , we 
use those of the most abundant ion, HCO + . Similarly, we regard the rate coefficients of Mg + as 
the representatives of those of metal ions M + . 

2.3.2. Reactions on Grain Surfaces 

Dust grains may gradually settle to a very thin layer about the midplane reducing the grain 
abundance in the major part of the disk (Nakagawa, Nakazawa, &; Hayashi 1981). We define a 
parameter f g that represents the abundance of dust grains relative to that in molecular clouds. 
The depletion parameter / g may be regarded as representing the evolutionary stage of the disk. 

For simplicity, we assume that all dust grains have the same size (radius) a, and that they 
consist of rocky and metallic materials. We take the inner density of grain material as 3 g cm -3 . 
We discuss the effects of the grain-size distribution and the effect of tiny grains in §3.3. We 
neglect icy material (H2, CH4, and NH3) in dust grains because they have only small effects on 
the ionization state of the disk (Umebayashi 1983). 

At the low temperatures concerned in this work, dust grains can adsorb ions and electrons. 
Thus, in addition to the radiative and dissociative recombinations of ions and free electrons in 
the gas phase, we have to take into account their recombination on grain surfaces, which occurs 
through grain-ion, grain-electron, and grain-grain collisions. As in the previous works (Nakano 
1984; Umebayashi & Nakano 1990; Nishi et al. 1991), we adopt the following model for the 
reactions on grain surfaces: 

1. When an ion hits a neutral grain or a positively charged grain, it sticks to the grain with a 
probability Si ~ 1. The sticking probability S e of an electron onto a neutral or negatively 
charged grain has been found to be between 0.2 and 1.0 (Nishi et al. 1991). We use S e = 0.6 
as in the previous works. 
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2. When an ion hits a negatively charged grain, or an electron hits a positively charged grain, 
the electron and the ion recombine and leave the grain surface using the energy released by 
the recombination. 

3. When grains of opposite charge collide in thermal Brownian motion, they neutralize 
themselves. We neglect coalescence of grains. 

We take account of the effect of electrical polarization of grains for the collision with an ion 
or an electron by using the fitting formula for the collision rate coefficient obtained by Draine & 
Sutin (1987). 



2.3.3. Ionization Processes 

In protoplanetary disks, charged particles are formed first by ionization of H2 and He by 
cosmic rays and radioactive elements. Table [3| lists the reaction rates of H2 and He with such 
energetic particles in units of the total ionization rate £ of a hydrogen molecule, which includes 
the ionization by secondary electrons. The ionization rate in the disk is effectively given by 



, x{r,z)\ ( £(r) - x(r,z) 

exp + exp 1 



+ Cr , (19) 



XCR / \ XCR 

where (cr ~ 1-0 x 10 'sec -1 is the ionization rate by cosmic rays in the interstellar space, and 
Cr, ~ 6.9 x 10 _23 / g sec -1 is the ionization rate by radioactive elements contained in the disk 
(Umebayashi & Nakano 1981). The depth 

POO 

X(r, z) = I p(r, z)dz (20) 

J z 

is the vertical column density from infinity to the position concerned. The attenuation length of 
the ionization rate by cosmic rays has a value xcr ~ 96 g cm~ 2 (Umebayashi & Nakano 1981). 
Because most of the radioactive elements are in rocky and metallic materials, the rate Cr is 
proportional to the depletion factor f g of dust grains. 



3. RESULTS 

3.1. The Minimum-Mass Solar Nebula — The Fiducial Model 

First, we investigate the stability of a disk rotating around a star of M* = 1 Mq with the 
surface density parameter /s = 1, which we presume to well represent the primitive solar nebula. 
We take p = 1/2 and To = 280 K in equation (|4|). We adopt the radius of dust grains a = 0.1 
/j,m and the grain depletion factor / g = 1 corresponding to early evolutionary stages where 
sedimentation of dust grains has hardly proceeded yet. We assume 61 = 0.2 and 82 = 0.02 for 
the parameters of element depletion, and we adopt fo 2 = 0.7 in accordance with the results of 
Mitchell, Kuntz, & Ginsburg (1978) on the calculation of molecular abundances. We call this the 
fiducial model. 
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3.1.1. Spatial Distribution of Charged Particles 

Figure [l] shows the distribution of relative abundances n(X)/nn of some representative 
particles X on the midplane for the fiducial model. The ionization rate at the midplane is also 
shown in the upper panel as a function of the radius r. Figure shows the resulting electrical 



conductivity a c given by equation (11). The abundances of charged particles and the electrical 
conductivity decrease as the distance from the central star decreases. When the time scale of 
magnetic field dissipation, 

H 2 

tdis ~ , (21) 

V 

is comparable to or smaller than the Keplerian orbital period £k = 2ir/Q, magnetic fields are 
almost decoupled from the gas. We call the distance from the central star at which tdis = £k holds 
the decoupling radius rd cc • Figure |3] depicts the ratio of the dissipation time tdis to the dynamical 
time tK as a function of r. We find tdis ^ *k at r < r^ cc 15 AU, where the hydrogen number 
density nn is higher than 10 12 cm -3 . 

In the region of r ^ rd ec > electrons and metal ions are the dominant charged particles as seen 
from Figure |]. Electrons and ions recombine mainly by collision with dust grains, and thus the 
relative abundances are inversely proportional to nn- The conductivity is mainly contributed by 
electrons, and is proportional to x(e). 

In the region of r ^ r^ cc , on the other hand, the number density of dust grains is higher than 
that of ions, and the dust grains with charge ±e are the dominant charged particles. Because the 
charge neutrality must be kept mainly by grains, or n(G(+l)) ~ n(G(— 1)), the number density 
of free electrons n(e) becomes considerably lower than that of metal ions n(M + ). The ratio 
n(e)/n(M + ) approaches the limiting value (m e /S e ™M) 1 '' 2 (Umebayashi 1983; Nakano 1984) as r 
decreases from 10 AU, where m e and myi are masses of electron and metallic ion, respectively. 
At r 7 AU, the column density of the disk exceeds the attenuation length of cosmic rays, XCR> 
and the abundance of charged particles reduces considerably due to decrease of the ionization 
rate. In the region of r <^ 2 AU, where the terrestrial planets orbit in the present solar system, 
the abundances of dust grains with electric charge ±e decrease as r decreases. This is due to 
recombination of adsorbed electrons and ions during grain-grain collision. At r <^ 1 AU, where the 
number density of ions is more than 10 4 times smaller than that of charged grains, the contribution 
to a c is dominated by charged grains as seen in Figure |. 

We next investigate the abundances of charged particles and the dissipation time of magnetic 
field off the midplane at some representative regions. Figure |] depicts the relative abundances 
of some representative particles as functions of the height \z\ from the midplane for the region 
r = 1 AU where the Earth may have formed. As \z\ decreases, the abundances of charged particles 
decrease, and as a result the electrical conductivity also decreases. At \z\ <^ 3.0 H, the number 
density of free electrons deviates from that of metal ions, and dust grains with electric charge dbe 
are the dominant charged particles. At \z\ ^ 1.0 H, the number density of charged grains decreases 
as \z\ decreases because of the charge neutralization by mutual collisions. The conductivity is 
mainly contributed by ions and electrons except at \z\ <^ 0.1 H where the number density of ions 
is lower than 10 -4 times that of charged grains. 
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Figure || shows the relative abundances of some representative particles for another typical 
region r = 30 AU where Neptune may have formed. Because of the low density nn <^ 10 10 cm~ 3 , 
electrons and metal ions are the dominant charged particles at any height. 

Figure |6| shows the dissipation time i^is as a function of \z\ at r = 1 and 30 AU. In the region 
r = 1 AU tdis is smaller than the dynamical time ix at \z\ ^ ^dec ~ 3.2 H where the number 
density nu exceeds 10 11 cm -3 . Notice that the column density of the well-coupled region, where 
^dis > % holds, is smaller than 10 -2 g cm -2 , which is much less than the attenuation length of 
cosmic rays XCR> and as a result ( as £cr as shown in the upper panel of Figure f|, different from 
Gammie (1996). The low ionization fraction in the decoupled regions is caused merely by the very 
high density, not by attenuation of cosmic rays. At r = 30 AU we find idis 3> at any z, and 
consequently the gas is strongly coupled with magnetic fields. 



3.1.2. Unstable Regions 

We shall find out the regions where the condition given by equation (||) is satisfied using the 
spatial distribution of magnetic diffusivity on the r-z plane obtained above. In such regions, the 
magnetorotational instability is operative (Paper I), and thus we call them the unstable regions. 

Since the stability condition also depends on the Alfven speed, we have to specify the spatial 
distribution of magnetic field strength. For the sake of simplicity, we consider only the vertical 
component of magnetic fields, and assume that the plasma beta at the midplane, f5 c = 2c^/v Ac , 
is constant in the disk, where vac is the Alfven speed at the midplane. Then, the distribution of 
magnetic field strength for the disks with q = 3/2, M* = 1 Mq, p = 1/2, and To = 280 K is given 
by 

Figures ^a and [7|b show the unstable regions for the cases of (3 C = 100 and 1000, respectively. 
In these figures the curve labeled Xid ea ,i/H = 1 represents the locus of the points at which the 
characteristic unstable wavelength in the ideal MHD limit, Ajdeals is equal to the scale height of 
the disk, 

\ idcal (r,z)^2n^^=H(r). (23) 

The curve with X res /H = 1 is the locus of the points at which the characteristic unstable 
wavelength in the resistive limit, A res , is equal to the scale height, 

\ ICS (r, z )^27T^^-=H(r). (24) 
v A (r,z) 

The critical condition for decoupling of magnetic field, idis = mentioned above is equivalent to 
equation when j3 c = 1. 



For a given radius, Ajdeai increases exponentially as the height \z\ increases because the Alfven 
speed becomes greater, while A res decreases because the magnetic diffusivity decreases as \z\ 
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increases. Therefore, the instability condition max(Aideai> X res ) <^ H [eq. @] is satisfied only in the 
striped area between the two critical curves Xidea\/H = 1 and X res /H = 1 in Figure [7|. 

The condition \\dc&\/H = 1 is easily reduced to 

8 \ 1/2 

which gives \z\ = 0.96 H and 1.8 H for B c = 100 and 1000, respectively. The magnetorotational 
instability is suppressed due to magnetic dissipation in the region below the critical curve 
X res /H = 1, which corresponds to Gammie's (1996) "dead zone". The dissipation process is less 
important for stronger magnetic field. Thus the stable region shrinks as the plasma beta decreases. 
At the points where X res /H = 1 holds, the hydrogen number density ran takes 0.8 - 1.7 x 10 11 
and 0.5 - 1.2 x 10 11 cm -3 for the cases of 8 C = 100 and 1000, respectively. At temperatures 
below several hundred Kelvin where thermal ionization is negligible, the magnetic diffusivity is 
determined mainly by the density and the ionization rate. Because the column density above the 
dead zone is far below the attenuation length of cosmic rays, the boundary of the dead zone is 
determined mostly by the local density in this case. 

The electrical conductivity is mainly contributed by electrons at the boundary of the dead 
zone. The relative abundance of electrons x(e) at the boundary is 1.1 - 2.7 x 10 -14 and 2.8 - 
7.1 x 10~ 14 for the cases of 3 C = 100 and 1000, respectively. The magnetic diffusivity is inversely 
proportional to x(e) in such a situation, and then we approximately have x(e) oc on the curve 
X Tes /H = 1 at a given radius. 

Let the critical radius r cr i t be the outer radius of the dead zone on the midplane, which is 
about 19 and 22 AU for the cases of 3 C = 100 and 1000, respectively. The region r ^ r cr i t is 
magnetorotationally unstable except at tenuous outer layers. In this region, because the magnetic 
field is amplified due to growth of fluctuations, the angular momentum can be transported 
effectively. As seen from Figure |7|b, the unstable region exists even at r < r cr i t above the dead 
zone. However, the column density of this unstable layer is so small that mass accretion cannot 
be as active as in the region r > r cr i t . Therefore, the accreting matter from the outer region must 
accumulate near the outer boundary of the dead zone, unless some other mechanisms of angular 
momentum transport are effective. 

In the following subsections, we shall depict the unstable regions on the r-z plane for various 



disk models. The criterion given by equation (24) is more important than the criterion (23) 



because the critical radius on the midplane r cr it is determined by equation (^J). Therefore, in the 
following we focus on the dead zone given by the criterion X TCS /H = 1. Because the results are not 
very sensitive to the plasma beta as seen from Figure |7], we show only the results for the cases of 
B c = 100 and 1000 in the following. 



3.2. Disks with Various Surface Densities 
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3.2.1. Dependence on the Disk Mass 

We investigate the disks with different surface densities rotating around a star of M* = 1 Mq. 
We consider the following cases of the surface density parameter; (a) /s = 0.3, (b) /s = 1 (the 
minimum-mass solar nebula), (c) /s = 3, and (d) /s = 10. Figure || shows the results for the cases 
with the depletion factor of dust grains / g = 1 corresponding to early evolutionary stages of the 
disks. The other parameters are the same as in the fiducial model. The solid and dashed curves in 
each panel represent the loci of X rcs /H = 1 for the cases of (5 C = 100 and 1000, respectively. Inside 
these curves are the dead zones where the ohmic dissipation suppresses the magnetorotational 
instability. The dotted curve shows the scale height of the disk z = H(r). 

Because the dead zone has a height close to the scale height H for most part of r < r crit , 
most of the matter at r < r cr i t is in the region where the instability is suppressed. Because the 
ohmic dissipation is not so efficient as to suppress the magnetorotational instability outside the 
dead zones in Figure ||, magnetic fields are amplified to some extent by the instability and angular 
momentum transport due to magnetic stress is expected though mass accretion through the layers 
above the dead zones may not be efficient because of their low column densities. 

We find that a higher mass disk is more stable for the magnetorotational instability; as 
the surface density parameter /s increases, the dead zone expands and the critical radius v^crit 
increases. This is because the recombination processes are more effective in denser regions, and 
thus the ionization fraction is lower in higher mass disks. In the case of (3 C = 100 (1000), the 
critical radius on the midplane, r cr ;t, is 13 (16), 19 (22), 26 (31), and 39 (45) AU for fc = 0.3, 1, 
3, and 10, respectively. At the points where \ ves /H = 1 holds, the hydrogen number density takes 
a value between a few x 10 10 and a few x 10 11 cm" 3 almost independent of fe. Therefore, the 
boundary of the dead zone is determined mainly by the local density ran; cosmic rays are hardly 
attenuated at these points as in the fiducial model. 

3.2.2. Dependence on the Power Index o/S 

So far we have considered only the disks with the power index of S, q = 3/2. There are no 
definite observations that constrain surface density profiles of protoplanetary disks. Some of the 
extrasolar planet candidates discovered by radial velocity oscillation orbit within just 0.3 AU, 
suggesting that they have migrated inward after their birth (e.g., Marcy, Cochran, and Mayor 
2000). Thus the mass distribution of the present-day solar system might be different from the 
primordial distribution. Here we examine some other models with different power index; (a) q = 2, 
(b) q = 3/2 (fiducial model), (c) q = 1, and (d) q = 0. We assume that all these models have equal 
mass 0.024 M Q between 0.1 and 100 AU. Thus the surface density at 1 AU is 4.8 x 10 3 , 1.7 x 10 3 , 
3.3 x 10 2 , and 6.6 g cm -2 for q = 2, 3/2, 1, and 0, respectively. The results for a disk of / g = 1 
around a star of M* = 1 Mq are shown in Figure ||. The other parameters are the same as in the 
fiducial model. The results are quite similar among these models though the critical radius ^crit 
takes somewhat different values. The critical radius on the midplane is mainly determined by the 
density which is about 1.6 x 10 11 and 1.1 x 10 11 cm -3 for j3 c = 100 and 1000, respectively. 
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3.3. Disks in Some Evolutionary Stages 

3.3.1. Depletion of Dust Grains 

The ionization fraction and the species of dominant charged particles are sensitive to the 
depletion factor of dust grains, f g , the abundance of dust grains relative to that in molecular 
clouds. As the dust grains in a protoplanetary disk sink gradually toward the midplane, f g 
decreases and the recombination rate on grain surfaces reduces except in a very thin layer around 
the midplane. As a result, the ionization fraction increases (Umebayashi & Nakano 1988) and the 
dead zone shrinks in the late evolutionary stages. 

Figure [l(] shows the relative abundances of some representative particles as functions of the 
height from the midplane at r = 1 AU for the case of f g = 1CP 4 . Electrons and metal ions are the 
dominant charged particles at any height z due to the low abundance of dust grains. Magnetic 
fields are well coupled with the gas in the larger part of the disk compared with the case of f g = 1. 
The electrical conductivity a c is mainly contributed by electrons, and is proportional to x(e) at 
any z. At \z\ <^ 1.0 H, the ionization by cosmic rays declines (upper panel of Fig. |l^) because the 
depth x(r, z) exceeds the attenuation length of cosmic rays, XCR- As seen from the dot-dashed 
curve in Figure ||, the ratio idisAK becomes less than unity at \z\ ^ Zd ec ~ 0.8 H, so that the gas 
decouples from magnetic fields in this region. This decoupled layer is significantly thinner than for 
the case of f g = 1. 

To clarify the effect of grain depletion, we investigate several cases with respect to the 
depletion factor of dust grains, (a) / g = 1 (fiducial model), (b) f g = HT 1 , (c) f g = 10" 2 , and (d) 
fg = 10 -4 , each of which might correspond to some evolutionary stage of the protoplanetary disk. 
Figure [ll| shows the results for the disk with the surface density parameter /s = 1 around a star 
of M* = 1 Mq. The other parameters are the same as in the fiducial model. 

We find that decrease of dust abundance makes smaller the region where magnetorotational 
instability is suppressed by ohmic dissipation. At later stages where only a very small fraction of 
the initial dust grains floats, magnetic fields are nearly frozen to the gas in most parts of the disk. 
This implies that the magnetorotational instability is operative at a larger part of the disk and 
the angular momentum transport is more effective in the later evolutionary stages. For the case of 
j3 c = 100, the critical radius r cr it is 19, 7.6, 3.8, and 2.0 AU for the models of f g = 1, HT 1 , 10~ 2 , 
and 10 -4 , respectively. 

For f g = 10~ 4 , the critical height z CT \t where X Tes /H = 1 holds is 1.0 H and 1.2 H for the cases 
of 0c = 100 and 1000, respectively, at 1 AU. At these points, electrons and metal ions are still 
dominant charged particles (see Fig. |l^). At \z\ = z cr jt in the region of r ~ 1 AU, the ionization 
rate is smaller than £cr by an order of magnitude. This is because the column density above the 
dead zone is nearly equal to XCR ~ 96 g cm -2 for these cases. 
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3.3.2. Growth of Grains 

The interstellar grains have a wide size distribution, e.g., between 0.005 and 0.25 /im (Mathis, 
Rumpl, & Nordsieck 1977; hereafter MRN). By coalescence of dust grains in protoplanetary disks, 
the size distribution and the relative number density of grains may change considerably. To clarify 
the effect of grain size on the unstable regions, we investigate the following four cases assuming 
for simplicity that all grains have the same radius; (a) a = 0.03 /im, (b) a = 0.1 /im (fiducial 
model), (c) a = 0.3 /im, and (d) a = 1 ftm. The relative number density of dust grains n g /nn is 
proportional to a -3 and is, for example, 3.3 x 10~ n and 8.9 x 10~ 13 for a = 0.03 and 0.1 /im, 
respectively. The results for a disk of /s = 1 and / g = 1 around a star of M* = 1 M Q are shown 



in Figure 12. The other parameters are the same as in the fiducial model. 



We find that the dead zone shrinks as the grain size increases. As the grains grow, the ion 
density increases because the total surface area of grains decreases and thus the recombination 
rate on grain surface decreases. For the case of (3 C = 100, the critical radius r cr it is 37, 19, 7.6, and 
4.6 AU for the models with a = 0.03, 0.1, 0.3, and 1 /im, respectively. 

For the case of a = 0.03 /im, the number density of dust grains becomes comparable to or 
higher than that of metal ions at nn ^> 10 cm . At the critical radius ^crit « 37 AU on the 
midplane, the abundance of metal ions x(M + ) ~ 1.4 x 10~ 12 is comparable to those of charged 
grains with ±e, though the conductivity a c is still mainly contributed by electrons with abundance 
x{e) ss 0.9 x 10~ 14 . 

So far we have assumed that all grains have the same radius. However, interstellar grains 
have a wide size distribution such as the MRN distribution for graphite and silicate grains, 

^^ = An R a~ 3 - 5 , a min < a < a max , (26) 
da 

where a min ss 0.005 /im, a max ss 0.25 /im, and A«1.5x 10~ 25 cm 2 - 5 (Draine & Lee 1984). For 
comparison, we show the results for this size distribution. Figure [l^ shows the vertical distribution 
of the relative abundances n(X)/nn of some representative particles at 1 AU. 

The qualitative features are essentially the same as those for the fiducial model (see Fig. ||). 
The total surface area of grains is contributed mainly by grains with radius a ~ a m i n in the MRN 
distribution. Therefore, decrease of a m j n enhances the recombination rate of ions and electrons 
on grains; a m i n can be smaller than 0.005 /jm (MRN). Dust grains with electric charge ±e are 
the dominant charged particles at the midplane everywhere within r ss 100 AU. In this case the 
magnetic field dissipation is more effective than in the fiducial model. The magnetorotational 
instability is suppressed in a wider region, ctncL the critical radii on the midplane rcrit 

= 44 and 

59 AU for the cases of (5 C = 100 and 1000, respectively, are 2 to 3 times larger than those for the 
fiducial model. 



4. DISCUSSION 
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4.1. Layered Accretion 



In paper I we have found that, in some disk conditions, the unstable modes have large 
amplitudes localized at the upper layer of the disk. This implies that the magnetorotational 
instability in the nonlinear regime enables the angular momentum transport only in this layer, 
and thus the layered accretion occurs as proposed by Gammie (1996). This layer corresponds to 
the unstable region at r ^ r crit shown by the stripes in Figure [7|. Since this unstable layer locates 
at a height a few times the scale height of the disk, the column density of this layer is very small. 
Therefore, the layered accretion might hardly affect dynamical evolution of the disks. However, 
the thickness of this layer is determined by the distribution of the magnetic diffusivity, which 
depends not only on the density distribution but also on the size and abundance of grains. 

First, we examine the dependence of the column density of the unstable layer, S uns , on the 
depletion factor of dust grains, / g . Because the column density of the layer above the critical 
curve \ide&\/H = 1 is quite small, we approximate S un s as the total column density of the disk 
subtracted by the column density of the dead zone. Figure |H| a shows E 

uns ei function of r for 

j3 = 100 for three cases of grain depletion, / g = 1, 10" 2 , and 10" 4 . The grain size is taken to be 
a = 0.1 ^m and the other parameters are the same as in the fiducial model. The column density 
S uns increases as the grain fraction decreases. Gammie (1996) assumed that the column density of 
the unstable region was comparable to the attenuation length xcr ~ 96 g cm~ 2 . We have found 
that the column density S uns is much smaller than xcr unless / g ^ 10 -4 . Magnetic fields can be 
decoupled from the gas even when cosmic rays are not attenuated. 

Next we examine the dependence of S uns on the grain size. Figure |l4| b shows the column 
density S uns for the models with a = 0.1, 0.3, and 1 /mi. Here we have assumed the grain fraction 
/ g = 1 and P c = 100. The column density S uns increases as the grain radius increases. However, 
S uns is at most a few g cm -2 within r ~ 1 AU, even when the grain is as large as a = 1 /j,m. 
Therefore, unless the abundance of dust grains is reduced to 10 -4 times the interstellar abundance, 
mass accretion through the unstable layer above the dead zone is very faint and the matter must 
accumulate around the outer boundary of the dead zone. 



We estimate the mass accumulation rate into the dead zone. Sano, Inutsuka, & Miyama 
(1998) investigated the nonlinear evolution of magnetorotational instability including the effect 
of ohmic dissipation using two-dimensional MHD simulations. Their results show that, at the 
saturated turbulent state, the efficiency of angular momentum transport depends on the strength 
of the initial magnetic field perpendicular to the disk, and that when the magnetic Reynolds 
number R m = v\/r]Q, is around unity, the a parameter of the viscosity (Shakura & Sunyaev 1973) 



4.2. 



Mass Accretion Rate 



is given by 




(27) 



where (5 is the plasma beta for the initial field. 
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Because the magnetic Reynolds number R m is close to unity around the critical radius, the 
mass accretion rate can be estimated with equation fl27|). The radial accretion velocity v r is given 
by v r ~ v/r, where v = ac s H is the kinematic viscosity. Then the mass accretion rate at v ~ t ct h 
is given by 

M ~ 2vrr crit S Wr « 4.9 x 10~ 7 /k (t^Xu) ^ (l^) ' Mq ^ ' (2§) 

for the disks whose structure is given by equations (||), (j5|), and (0) with q = 3/2 and M* = 1 Mq. 
For the case of (5 C = 1000, the mass accretion rate at r ~ r cr i t ~ 22 AU is an order of magnitude 
smaller. 

The typical age of T Tauri stars is estimated to be 10 6 yr (Kenyon & Hartmann 1995), and 
the lifetime of protoplanetary disks must be comparable to or longer than the age of T Tauri 
stars. If mass accretion with a rate M ~ 10 -6 Mq yr -1 continues and the accreted matter is 
distributed relatively widely inside r ~ r cr it, the disk of this part becomes gravitationally unstable 
at some stage, and the resultant gravitational torque causes mass accretion in this region. This 
makes the accretion onto the central star non-steady. If the accreting matter accumulates near the 
critical radius r ~ r cr i t and a dense ring forms, a planet or a brown dwarf would be formed due to 
fragmentation of the ring (Nakano 1991). In this case, the critical radius gives the position where 
the first small companion forms. 

Higher mass disks have higher accretion rate as seen from equation (|28|). The mass accretion 
rate for a disk of /s = 10 at r ~ r cr - lt ~ 39 AU is 7 times larger than that of the fiducial model at 
r cr it ~ 19 AU. For the model of /s = 0.3, we have the mass accretion rate M ~ 1.8 x 10~ 7 Mq yr -1 
at r ~ r cr it ~ 13 AU. Therefore, higher mass disks become gravitationally unstable in shorter time 
scale. 



4.3. Evolution of Dust Grains in Protoplanetary Disks 

Investigating the growth and sedimentation of dust grains in protoplanetary disks, Nakagawa 
et al. (1981) showed that dust grains can grow by sticking through mutual collisions in thermal 
Brownian motion, e.g., when the grain size is smaller than microns at r ~ 1 AU. The collision time 
in Brownian motion, t co \, is given by 

, _ 2 (P S /3) 3/2 a 5/2 _ 21 / P S \ 3/2 ( a \ 5/ V C g V7 r Y /2 f f29l 
tcol ~ 27r Q P {kT) 1 ' 2 ~ 21 Ul^J \U^) {OMU) [iMj) tK ' (29) 

where p g and a are the internal density and the mean radius, respectively, of dust grains, and Cg is 
the mass fraction of dust grains, and the last expression is at the midplane of the fiducial model. 
As seen from equation (^9|), dust grains in dense regions grow to micron size in a time comparable 
to the dynamical time scale, t^, and much shorter than the time scale of sedimentation, 

W = tK = 5.2 x 10 5 ( 1 ( -i-) _1 yr , (30) 

Z7r p g a \6 g cm J \1 pm / 

where the last expression is at the midplane of the fiducial model and is independent of r. 
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Numerical calculations of Nakagawa et al. (1981) show that the mass fraction of floating 
grains, / g , at r = 1 AU is ~ 1CT 1 , 10~ 2 , and 10~ 4 at the time t~ 2 x 10 3 , 10 4 , and 10 5 yr, 
respectively, and that the size of floating grains, a, is 1 - 10 fim in the late stages. Although 
laminar flow was assumed in their analysis, we discuss the evolution of the unstable regions using 



their results. Figure 12 d shows the unstable region for the model with a = 1 //m, an order of 
magnitude larger than in the fiducial model. The critical radius on the midplane r cr it is 4.6 and 5.8 
AU for the cases of (3 C = 100 and 1000, respectively. The dead zone shrinks as the sedimentation of 
dust grains proceeds; in the cases of / g = 10 _1 , f g = 10 -2 , and / g = 10 -4 , the critical radius of the 
models with a = 1 is 3.2 (3.8), 2.3 (2.7), and 1.4 (1.6) AU, respectively, for /3 C = 100 (1000). 
In all these cases the critical radius r CT [ t is less than 6 AU. If dust grains in protoplanetary disks 
could grow to micron size, most parts of the disk become magnetorotationally unstable. Therefore, 
sedimentation of dust grains must be prevented by turbulence induced by this instability contrary 
to the assumption of Nakagawa et al. (1981). 

According to the standard scenario of solar system formation (Hayashi et al. 1985), planets 
form by accretion of planetesimals, which are produced through gravitational instability of a thin 
dust layer (Goldreich & Ward 1973). The thin layer is formed by sedimentation of dust, which 
must occur only in a quiescent disk (Mizuno, Markiewicz, & Volk 1988). By sedimentation the 
disk becomes magnetorotationally unstable and turbulent, preventing sedimentation as shown in 
§3.3. Therefore we have to investigate the gas motion and the evolution of dust simultaneously. 



4.4. Ionization by X-Rays 

So far we have considered only cosmic rays and radioactive elements as the ionization 
sources of protoplanetary disks. In this subsection we consider the effect of photo-processes. 
The protoplanetary disks are exposed to the UV radiation from the central star and from the 
interstellar space. Because of the very small attenuation length, a few x 10~ 3 g cm~ 2 , the 
ionization by UV radiation is inefficient except for very thin surface layers. 

Young stellar objects are strong X-ray sources (e.g., Glassgold, Feigelson, & Montmerle 2000). 
The X-ray luminosity of low-mass young stellar objects is typically in the range Lx ~ 10 28 - 10 30 
erg s , or a factor of 10 2 - 10 3 above the contemporary solar levels. Some X-ray sources exhibit 
high-amplitude rapid flares with peak luminosity Lx ~ 10 30 - 10 32 erg s . These emissions 
require the presence of a large volume of high-density plasma at a temperature of about 10 K. 
The plasma must be magnetically confined, probably in large loops on a scale comparable to 
or larger than the X-ray emitting star itself. Field lines linking the star with the disk at the 
corotating radius may explain the strongest of these powerful flares (Shu et al. 1994; Hayashi, 
Shibata, & Matsumoto 1996). The X-rays can induce a wide variety of changes in chemical and 
physical properties of protoplanetary disks. 

Soft X-rays can be the most promising ionization sources for protoplanetary disks as well as 
cosmic rays. While the underlying ionization mechanism is due to electronic collisions in both 
cases, X-rays are absorbed in a smaller column of matter than cosmic rays. Thus cosmic rays 
provide a more global, low-level ionization of the disks, whereas X-rays can produce a localized 
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high-level ionization. 

Igea & Glassgold (1999) investigated the X-ray transfer and ionization in the disks around 
young stellar objects using a Monte Carlo method. The X-ray emitting region was modeled as 
a ring of radius r = 10 Rq at height z = 10 R@ with Lx = 10 29 erg s _1 based on the x-wind 
model of Shu et al. (1994). They obtained the ionization rates at 1 - 10 AU of the minimum-mass 
solar nebula as functions of the vertical column density into the disk by hydrogen number, 
N±. Their results show that the ionization rate by X-rays is larger than that by cosmic rays, 
Ccr ~ 10~ 17 s _1 , only at N± smaller than 10 25 to 2 x 10 23 cm -2 depending on r, which is much 
smaller than the total column density of the disk at each r. Therefore, most of our results are 
preserved even if X-ray ionization is taken into account. For instance, the thickness of the dead 
zone is hardly affected by the X-ray ionization. Even if cosmic rays are excluded from the inner 
regions of protoplanetary disks by, e.g., stellar winds (Parker 1960), the active surface layer can be 
maintained by the X-ray ionization (Igea & Glassgold 1999). Igea k, Glassgold (1999) neglected 
the effects of dust grains on ion densities and on interaction with magnetic fields. However, the 
column densities of unstable regions Suns must be affected by the characters of dust grains in the 
disk in the same way as in the cases with the cosmic ray ionization investigated in §3. 



4.5. Uncertainties in the Disk Temperature 

The temperature distribution given by equation (||) with p = 1/2 holds for the disk transparent 
to the radiation from the central star. Actually, the inner dense regions of protoplanetary disks 
would be optically thick and there would be some other heating sources such as viscous heating 
and dissipation of magnetic fields. However, the results are not very sensitive to the temperature 



as long as it is below several hundred Kelvin where the thermal ionization is inefficient. Figure |15 
shows the column densities of the unstable layers for the models with temperature ten times higher 
(dashed curve) and ten times lower (dot-dashed curve) than the fiducial model (solid curve). The 
other parameters are the same as those of the fiducial model. The difference in E uns is surprisingly 
small for such large changes in the temperature. We also investigated for p = and 3/4, and 
obtained quite similar results assuming as above that the thermal ionization is inefficient. When 
the temperature is higher than several hundred Kelvin, the thermal ionization becomes efficient 
and there would be significant effects on the instability. This would happen only in the innermost 
region r ^ 0.1 AU for the minimum-mass model, which is much smaller than the dead zones 
obtained in this paper. Notice that the temperature of steady accretion disks (p = 3/4) with 
typical mass accretion rate is within the range we examined here. 



4.6. Validity of Ionization-Recombination Equilibrium 

We have used the densities of charged particles obtained by assuming that the gas is in 
ionization-recombination equilibrium. This assumption is valid when the relaxation time to 
an equilibrium state for each kind of charged particles is much shorter than the characteristic 
time scale of protoplanetary disks. The relaxation time i r (Xj) to the ionization-recombination 
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equilibrium for a species Xj is approximately given by the minimum of 

n(Xi) 



/3 ijk n(Xj)n(X k 

We estimate i r (Xj) by using the equilibrium abundances of particles adopted in the previous 



(31) 



section. Figure |l6| a shows the time scales near the midplane as functions of r, and Figure 16b 
shows those in the Earth's region (r = 1 AU) as functions of \z\, for the fiducial model. 



As seen from Figure |l6p, the gas near the midplane is almost in ionization-recombination 
equilibrium at r ^ 1.3 AU because t r is smaller than tx for all the particles. At r <J 1 AU, 
however, t T for neutral grains is much longer than tx- Because t T for ions and electrons is much 
shorter than t r for grains, ions and electrons always take quasi-equilibrium abundances determined 
for the temporary charge-state distribution of grains. When the charge-state distribution deviates 
significantly from the equilibrium, the relaxation time to its true equilibrium is about 250 yr, much 
longer than tx (see Fig [i~6|a) . However, because this is much shorter than the sedimentation time 
of grains in this region of the solar nebula (Nakagawa et al. 1981), the ionization-recombination 
equilibrium can be attained even in this region. 

In Earth's region the time scale t r for metal ions and electrons is longer than tx only at 
\z\ ;> 3.7 H as seen from Figure |l6|b. Because most of the disk matter is inside this layer, a possible 
deviation from the ionization-recombination equilibrium hardly affects the results obtained in the 
previous sections and subsections. 

As the depletion factor f g of dust grains decreases, the time scale t r for dust grains decreases 
considerably in regions where ions and electrons are the dominant charged particles. For ions 
and electrons t r remains at least 10 times smaller than tx even at such late evolutionary stages 
as fg = 10 -4 . In the inner regions where dust grains are the dominant charged particles, t r for 
dust grains increases extensively as r decreases, and it finally exceeds tx at a certain position. For 
the case of / g = 10 -4 the time scale t r for charged grains exceeds tx only at r ^ 0.66 AU, and 
its value is at most 300 times longer than tx even in the innermost region r <J 0.5 AU. Because 
the evolutionary time scale of the protoplanetary disk is much longer than tx, we can regard that 
the solar nebula is in ionization-recombination equilibrium even at stages of / g ~ 10~ 4 except 
at high \z\ where little matter exists. The situation is qualitatively the same for the other disks 
investigated in the previous section. 



5. SUMMARY 

We have investigated the ionization state and magnetorotational instability for various 
models of protoplanetary disks taking the effect of charged grains into account. Our findings are 
summarized as follows. 

1. In early evolutionary stages of the minimum- mass solar nebula, we have found that 
magnetorotational instability operates outside the critical radius T cr n « 20 AU. Inside the 
critical radius most part of the disk is stabilized by ohmic dissipation. Thus, the accreting 
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matter accumulates near the critical radius unless some other mechanisms of angular 
momentum transport are effective. 

2. Just outside this critical radius, the mass accretion rate is 10 -7 - 10 -6 M Q yr -1 . By- 
accumulation of matter at this rate the inner region becomes gravitationally unstable within 
the typical lifetime of protoplanetary disks. This suggests that non-steady accretion onto 
the central star occurs as in FU Orionis objects. 

3. The critical radius r crit depends on the mass of the disk. Because the abundances of charged 
particles are lower in denser regions, the disks of higher mass have larger stable regions and 
thus larger r cr i t . 

4. As dust grains in the disk grow or sediment toward the equatorial plane, the critical 
radius r cr it decreases. This suggests that magnetorotational instability is important in late 
evolutionary stages of the disks exciting turbulence and suppressing sedimentation of dust 
grains. Therefore, to clarify the processes of planet formation, behavior of grains should be 
studied in connection with gas motion. 

5. Even in the inner regions r <J r cr i t , some part of the surface layer is magnetorotationally 
unstable. However, because the column density of this layer is so low unless the grain 
fraction is as low as f g <; 10~ 4 that the layered accretion must not be effective at least in the 
early evolutionary stages. 



Numerical computations were partly carried out at the Astronomical Data Analysis Center of 
the National Astronomical Observatory, Japan. This work was partly supported by Grant-in-Aid 
for Scientific Research on Priority Areas (10147105) of the Ministry of Education, Science, Sports, 
and Culture of Japan. 



MOMENTUM TRANSFER RATE BETWEEN ELECTRONS AND 

NEUTRAL PARTICLES 



M. Hayashi (1981) compiled the values of the momentum transfer cross sections of an electron 
in collision with a hydrogen molecule and a helium atom at low energies, which are not much 
different from the geometrical cross sections and are significantly smaller than the Langevin's 
values. From these values we have obtained the empirical formulae for the momentum transfer 
rate coefficients as functions of the temperature T as follows, 



{av) 



(av) 



e-H 2 



c-Hc 



irm J 
fSkTX 1 / 2 



7.94 x 10~ 16 - 2.52 x 10~ 16 logT + 1.19 x 10~ 16 (logT) s 



V vrm 

where m is the reduced mass. 



5.30 x 10~ 16 -4.70 x 10~ 16 log T + 2.31 x 10" 16 (log Tf 



(Al) 



(A2) 
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Equations (Al) and ( |A2| ) give the rate coefficients about two orders of magnitude smaller 
than those used previously by many authors. Thus, we find that electrons are much more tightly 
coupled with magnetic field than previously considered. 
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Fig. 1. — Relative abundances n(X)/nn (lower panel) of some representative particles on the 
midplane of a disk with q = 3/2, / E = 1, p = 1/2, T = 280 K, M* = 1 M , and f g = 1 as 
functions of the distance r from the star. We call this the fiducial model because it is expected to 
well represent the primitive solar nebula in an early evolutionary stage. We adopted the depletion 
factors of heavy elements from the gas phase 5\ = 0.2 and 82 = 0.02, the fraction of gas phase 
oxygen in the form of O2, H2O, or OH /o 2 = 0.7, and the radius of dust grains a = 0.1 //m. 
G(n) represents grains with electric charge ne. The upper panel shows the ionization rate £ on the 
midplane as a function of r. 
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Fig. 2. — The electrical conductivity a c on the midplane of the fiducial model (same as in Fig. |l|) 
as a function of the distance r from the star. Also shown are the contributions of electrons, ions, 
and charged grains to the conductivity. 
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Fig. 3. — The ohmic dissipation time of magnetic field, t^g, on the midplane relative to the 
Keplerian orbital period ix as a function of r for the fiducial model. The magnetic diffusivity r\ 
given by equation ([H]) is evaluated with the abundances of charged particles shown in Figure |l[ 
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Fig. 4. — Relative abundances n(X)/nn of some representative particles as functions of the height 
\z\ from the midplane at r = 1 AU for the same disk as in Figure |l|. The upper panel shows the 
ionization rate £ as a function of \z\. 
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Fig. 5. — Relative abundances n(X)/nn of some representative particles as functions of the height 
\z\ from the midplane at r = 30 AU for the same disk as in Figure |]. The upper panel shows the 
ionization rate £. 
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Fig. 6. — The ratio of the ohmic dissipation time idis to the Keplerian orbital period £k as a function 
of \z\ at r = 1 AU {solid curve) and 30 AU {dashed curve). The magnetic diffusivity is evaluated 
with the abundances of charged particles shown in Figures |I| and || for the fiducial model. The 
dot-dashed curve is for the case of f g = 10 -4 at r = 1 AU. 
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Fig. 7. — The unstable regions (striped) for the cases of (a) the plasma beta at the midplane 
j3 c = 100 and (b) (3 C = 1000 for the fiducial model (same as in Fig. @). The dotted curve denotes 
the scale height of the disk, z = H(r). In the region above the curve \ide&\/H = 1, the instability 
does not effectively work because the wavelengths of the unstable modes are larger than the scale 
height of the disk. The region below the curve X Tes /H = 1 is stable because the magnetic dissipation 
is effective, and corresponds to Gammie's (1996) "dead zone". 
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Fig. 8. — Unstable regions in the disks with various surface densities; (a) /s = 0.3, (b) /s = 1 (the 
minimum-mass solar nebula), (c) /s = 3, and (d) /s = 10. The other parameters are the same as 
those of the fiducial model. The solid and dashed curves represent the loci of X TCS /H = 1 for the 
cases of the field strength [3 C = 100 and 1000, respectively, inside of which is magnetorotationally 
stable. The dotted curve represents the scale height of the disk z = H(r). 
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Fig. 9. — Unstable regions in the disks with various power indices for the surface density; (a) q = 2, 
S = 4.8 x 10 3 g cm' 2 , (b) q = 3/2, S = 1.7 x 10 3 g cm" 2 , (c) q = 1, S = 3.3 x 10 2 g cm" 2 , 
and (d) q = 0, So = 6.6 g cm~ 2 . The other parameters are the same as those of the fiducial 
model. All these models have the same disk mass 0.024 M between 0.1 and 100 AU. The solid 
and dashed curves represent the loci of X rC s/H = 1 for the cases of the field strength (5 C = 100 and 
1000, respectively, inside of which is magnetorotationally stable. The dotted curve represents the 
scale height of the disk z = H(r). 
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Fig. 10. — Relative abundances n(X)/nn of some representative particles for the case of f g = 10~ 
as functions of the height \z\ from the midplane at r = 1 AU. The other parameters are the sam 
as those of the fiducial model. The upper panel shows the ionization rate (. 
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Fig. 11. — Unstable regions for the disk of various evolutionary stages; (a) f g = 1, (b) f g = 10 _1 , 
( c ) fg = 10 -2 , and (d) / g = 10 -4 . The other parameters are the same as those of the fiducial model. 
The solid and dashed curves represent the loci of X res /H = 1 for the cases of the field strength 
(3 C = 100 and 1000, respectively, inside of which is magnetorotationally stable. The dotted curve 
represents the scale height of the disk z = H(r). 
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Fig. 12. — Unstable regions for the disks with various grain sizes; (a) a = 0.03 /jm, (b) a = 0.1 
^m, (c) a = 0.3 /im, and (d) a = 1 //m. The other parameters are the same as those of the fiducial 
model. The solid and dashed curves represent the loci of X TCS /H = 1 for the cases of the field 
strength (3 C = 100 and 1000, respectively, inside of which is magnetorotationally stable. The dotted 
curve represents the scale height of the disk z = H(r). 
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Fig. 13. — Relative abundances n(X)/nn of some representative particles as functions of the height 
\z\ from the midplane at r = 1 AU for the MRN size distribution of grains. The other parameters 
are the same as those of the fiducial model. The upper panel shows the ionization rate (. 
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Fig. 14. — The column density of the unstable layer, S uns , (a) for the models of grain depletion 
fg = 1, 10 -2 , and 10 -4 , and (b) for the models of grain radius a = 0.1, 0.3, and 1 /im. The other 
parameters are the same as those of the fiducial model. 
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Fig. 15. — The column density of the unstable layer, S uns , for the models with temperature 10 
times higher (dashed curve) and 10 times lower (dot-dashed curve) than the fiducial model (solid 
curve) . The other parameters are the same as those of the fiducial model. 
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Fig. 16. — The relaxation time t r to the ionization-recombination equilibrium for some 
representative particles (a) on the midplane as functions of r, and (b) at r = 1 AU as functions of 
\z\, for the fiducial model. The Keplerian orbital period tn is also shown by the dotted line. 
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Table 1. Abundance of Elements and their Fraction in the Gas Phase 



Element 


Abundance 


Chemical Species 


Fraction in Gas Phase 


H 


1 


H 2 


1 


He 


9.75 x 1CT 2 


He 


1 


C 


3.62 x 1CT 4 


CO 


6i 





8.53 x 1CT 4 


o 2 ,o 


Si 


Metal Total 


7.97 x 1CT 5 


M 


s 2 



Table 2. Two-Body Reactions and their Rate Coefficients 









Reaction 








Rate coefficient (cm 3 s x ) 


H+ 


+ 


Mg 


-» Mg+ 


+ 


H 




1.10 x 10~ 9 


H+ 


+ 





- 0+ 


+ 


H 




7.00 x 10- 10 cxp(-232.0/T) 


H+ 


+ 


o 2 


- 0+ 


+ 


H 




2.00 x 10~ 9 


HJ 


+ 


H 2 


- H+ 


+ 


H 




2.08 x 10~ 9 




+ 





-» OH+ 


+ 


H 2 




8.00 x 10~ 10 


H 3 + 


+ 


Mg 


-> Mg+ 


+ 


H + 


H 2 


1.00 x 10~ 9 


H 3 + 


+ 


CO 


-» HCO+ 


+ 


H 2 




f .70 x 10~ 9 


H 3 + 


+ 


o 2 


-» 2 H+ 


+ 


H 2 




5.00 x 10- 9 exp(-150.0/T) 


Hc+ 


+ 


H 2 


-» H+ 


+ 


H + 


He 


3.70 x fO" 14 exp(-35.0/T) 


He+ 


+ 


CO 


- C+ 


+ 


+ 


He 


1.60 x 10~ 9 


Hc+ 


+ 


o 2 


- 0+ 


+ 


+ 


He 


1.00 x 10~ 9 


C+ 


+ 


Mg 


-» Mg+ 


+ 


c 




1.10 x 10~ 9 


C+ 


+ 


H 2 


- CH+ 


+ 


hv 




4.00 x 10- 16 (T/300)-°- 20 


c+ 


+ 


o 2 


-» C0+ 


+ 







3.80 x 10~ 10 


c+ 


+ 


o 2 


-> 0+ 


+ 


CO 




6.20 x 10~ 10 


HCO+ 


+ 


Mg 


-> Mg+ 


+ 


HCO 




2.90 x 10~ 9 


H+ 


+ 


e 


-» H 


+ 


hv 




3.50 x 10- 12 (T/300)-°- 75 


H 3 + 


+ 


e 


- H 2 


+ 


H 




1.50 x 10~ 8 (T/300)-°- 50 


H 3 + 


+ 


e 


-» H 


+ 


H + 


H 


1.50 x 10~ 8 (T/300)- - 50 


He+ 


+ 


e 


-» He 


+ 


hv 




2.36 x 10- 12 (T/300)~°- 64 


C+ 


+ 


e 


-» C 


+ 


hv 




4.40 x 10- 12 (T/300)-°- 61 


Mg+ 


+ 


e 


-» Mg 


+ 


hv 




2.80 x 10- 12 (T/300)~°- 86 


HCO+ 


+ 


e 


-» CO 




H 




1.10 x 10~ 7 (T/300)" 1 - 00 
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Table 3. Ionization Rates 



I 


leaction 






Rate 


H 2 


- H 2 + 


+ 


e 


0.97( a 


H 2 


-» H+ 


+ 


H + e 


0.03C 


He 


-» Hc+ 


+ 


e 


0.84C 



a £ is the total ionization rate for a hydrogen 
molecule. 



